Spectral Gaps for a Metropolis–hastings Algorithm in Infinite Dimensions By
نویسندگان
چکیده
We study the problem of sampling high and infinite dimensional target measures arising in applications such as conditioned diffusions and inverse problems. We focus on those that arise from approximating measures on Hilbert spaces defined via a density with respect to a Gaussian reference measure. We consider the Metropolis–Hastings algorithm that adds an accept– reject mechanism to a Markov chain proposal in order to make the chain reversible with respect to the target measure. We focus on cases where the proposal is either a Gaussian random walk (RWM) with covariance equal to that of the reference measure or an Ornstein–Uhlenbeck proposal (pCN) for which the reference measure is invariant. Previous results in terms of scaling and diffusion limits suggested that the pCN has a convergence rate that is independent of the dimension while the RWM method has undesirable dimension-dependent behaviour. We confirm this claim by exhibiting a dimension-independent Wasserstein spectral gap for pCN algorithm for a large class of target measures. In our setting this Wasserstein spectral gap implies an L2-spectral gap. We use both spectral gaps to show that the ergodic average satisfies a strong law of large numbers, the central limit theorem and nonasymptotic bounds on the mean square error, all dimension independent. In contrast we show that the spectral gap of the RWM algorithm applied to the reference measures degenerates as the dimension tends to infinity.
منابع مشابه
Approximating Bayes Estimates by Means of the Tierney Kadane, Importance Sampling and Metropolis-Hastings within Gibbs Methods in the Poisson-Exponential Distribution: A Comparative Study
Here, we work on the problem of point estimation of the parameters of the Poisson-exponential distribution through the Bayesian and maximum likelihood methods based on complete samples. The point Bayes estimates under the symmetric squared error loss (SEL) function are approximated using three methods, namely the Tierney Kadane approximation method, the importance sampling method and the Metrop...
متن کاملModified Metropolis-Hastings algorithm with Delayed Rejection
The development of an efficient MCMC strategy for sampling from complex distributions is a difficult task that needs to be solved for calculating small failure probabilities encountered in high-dimensional reliability analysis of engineering systems. Usually different variations of the Metropolis-Hastings algorithm (MH) are used. However, the standard MH algorithm does generally not work in hig...
متن کاملModified Metropolis-Hastings algorithm with Delayed Rejection for High-Dimensional Reliability Analysis
The development of an efficient MCMC strategy for sampling from complex distributions is a difficult task that needs to be solved for calculating small failure probabilities encountered in high-dimensional reliability analysis of engineering systems. Usually different variations of the Metropolis-Hastings algorithm (MH) are used. However, the standard MH algorithm does generally not work in hig...
متن کاملDiffusion Limits of the Random Walk Metropolis Algorithm in High Dimensions by Jonathan
Diffusion limits of MCMC methods in high dimensions provide a useful theoretical tool for studying computational complexity. In particular, they lead directly to precise estimates of the number of steps required to explore the target measure, in stationarity, as a function of the dimension of the state space. However, to date such results have mainly been proved for target measures with a produ...
متن کاملDiffusion Limits of the Random Walk Metropolis Algorithm in High Dimensions
Diffusion limits of MCMC methods in high dimensions provide a useful theoretical tool for studying computational complexity. In particular, they lead directly to precise estimates of the number of steps required to explore the target measure, in stationarity, as a function of the dimension of the state space. However, to date such results have mainly been proved for target measures with a produ...
متن کامل